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We compute the correlation functions of the eigenvalues in the Gaussian unitary ensemble using the 
fermionic replica method. We show that non-trivial saddle points, which break replica symmetry, 
must be included in the calculation in order to reproduce correctly the exact results for the corre- 
lation functions at large distance. 
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I. INTRODUCTION 

Random Matrix Theory (RMT) has found broad applications in physics ranging from nuclear spectra to electrons in 
metallic-grains, see Refs.Eln for reviews. The first mathematically rigorous results for the level statistics were derived by 
GaudinB and DysonQ using the method of orthogonal polynomials. The progress in the theory of Anderson localization 
in theJate seventies established a .close relation between the RMT and matrix cr-models initially in their replicated 
fornJj'B. Despite of a lot of effortsLiB, however, the initial attempts to reproduce Dyson's results from the replicated 
cr-models were not successful: while the density of eigenvalues could be found easily, it seemed that the eigenvalue 
correlations could not be obtained from the replicated cr-models. The most detailed account of such attempts was 
probably given by Zirnbauer and Verbaarschotu, who computed both bosonic and fermionic replicated cr-models and 
obtained different results, both-differing from the correct one. Only the supersymmetric (SUSY) formulation of the 
cr-model introduced by Efetov£3 gave a correct-and beautiful- way of calculating the correlations of eigenvalues from a 
cr model formulation. It has thus become a common wisdom for fifteen years that the SUSY is the only field theoretic 
method capable to compute the RMT level statistics, while the internal subtleties of the replica method seem to make 
it inapplicable for this task. 

Looking at the replicated sigma model approach, it is clear that there is one underlying assumption in the existing 
computations, which is the absence of any spontaneous breaking of replica symmetry. The study of the statistics 
of eigenvalues of a large N x N random matrix is mapped exactly onto a cr-model, where the action is of order 
N. Among various saddle points of the cr-model which could a-priori contribute, only the trivial 'replica symmetric' 
one was discussed. In this paper we revisit the problem, and consider all possible saddle points. Because of the 
symmetry of the cr-model the saddle points are actually saddle point manifolds. We show that the computation of 
the r-point correlation function involves 2 r saddle point manifolds. For the one-point function (the level density), 
the trivial saddle point gives Wigner's semi-circle law, whereas the second one contributes to order 1/N, and is, in 
fact, needed to obtain the oscillatory component of the density of states (DOS). For the two-point function, the effect 
is even more dramatic since one of the extra three non-trivial saddle point manifolds contributes to the leading order 
at large N (the other two being 1/N corrections). Taking it into account gives the correct result for -the two-point 
correlation function. A similar situation was actually found in the SUSY, where Andreev and Altshulertil showed that 
the asymptotic behavior of the two level correlations at large energy differences (in units of the mean level spacing) 
can be obtained by the saddle point evaluation of the SUSY cr-model, including one extra non-trivial saddle point. 

We restrict ourselves here to the Gaussian unitary ensemble (GUE), where we derive the exact results for the 
DOS, including oscillatory 1/N correction as well as the large energy difference behavior of the two-point correlation 
function using the fermionic replica formalism. The strategy is simple enough to be easily generalized for the higher 
correlation functions for which it also leads to known GUE expressions. Technically, a very useful step is to integrate 
exactly over the angular degrees of freedom of the cr-model, using the Itzykson-Zuber integrafia, which leaves the n 
eigenvalues of the replicated cr-model as the only integration variables. The saddle points are thus discussed on the 
level of eigenvalues, and the resulting replica symmetry breaking appears to be a particularly simple version of the 
vector replica symmetry breaking mechanism encountered in several disordered systemsEj. 

The paper is organized as follows: in section [n] we present the calculations of the average DOS, including the 
1/N corrections. This section is also used to introduce notations and illustrate the technique on this simple, but 
instructive example. Calculation of the two-point correlation function is given in section |ll]. Finally sections |y| and 
are devoted to discussions of the method and remaining open questions. Some technical details are presented in 
three appendices. 
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II. DENSITY OF STATES 



We are interested in the spectral properties of the random N x N Hermitian Hamiltonians, H, with the Gaussian 
probability distribution function 



P(H) = 2^" 1 )/ 2 ^ /2 ex p{-f Trff 2 } 
We begin with the calculation of the one-point function Si (E) defined as 



(1) 



Si{E) = N- lr Tr(E - H)- 1 , (2) 

where the complex energy E has an infinitesimal negative imaginary part. The bar stands for the averaging over H 
with the measure given by Eq. ([[J). The DOS is given by the imaginary part of this correlation function 

v(E) = -ir~ 1 QSi(E) . (3) 

Introducing 2N anti-commuting variables ip x ,ip x , where a; is a discrete index x G {1, . . .N}, the correlation function 
may be written as 



S 1 (E) = N- 1 -^hxZ{E) : (4) 

where 

r n ( n } 



Z{E)= J n#-#« ex p{- E ^AES xy -H xy )^ y \ 

x—1 \ x^y—1 ) 



We use now the replica trick to write the logarithm as InZ = lim„^o(Z" — l)/n. As a result the correlation function 
takes the form 

1 o(™)/tn . o(")/cn _ at- i ^ 



Sr(E) = Um -Sr(^); S^ 1 {E) = — Z(Ey , (6) 

n^O n oE 

It is convenient to introduce the generating function Z^ n \E) = Z(E\) . . . Z(E n ), where E is a diagonal n x n matrix 
which has the form E — diag{£i, ... ,E n }. In the limit E\,... E n — > E the generating function Z^ n \E) approaches 
Z(E) n and gives the one-point function according to Eq. (^). The generating function Z^ n \E) may be written as the 

fermionic integral analogous to that of Eq. (|5|). The integration runs now over 2nN fermionic variables ip x , ijf x , where 
x £ {1 . . . N} and j £ {1 . . . n}. Performing the averaging over the random matrix H one finds 



- N n j ^ n I 

/ II II 0**4 exp I -YjlE^l - — T jk T kj \ 

x=lj=l ^ x,j j-k=l J 



ZW(E) = I Mil'/' •;' <U. -j. <-xp { -) j ' J i E jI j , - — > ^ T lh T, , } . (7) 

where Tj k = ^2 x ip x tp x . Introducing the Hubbard-Stratonovich decoupling of the last term with the x-independent, 
Hermitian n x n matrix Q and integrating out the anti-commuting variables, one obtains the er-model in the formEi: 

ZM(E) = c n J d[Q] exp | - y Tr Q 2 + N Tr ln(£ - i Q) j , (8) 



C «=(£) 2"(- 1 )/ 2 . (9) 

The standard route to study the large N limit of the DOS prescribes to take all energies equal to E, look for the 
saddle points of the functional integral, Eq. (||), and then consider fluctuations around the saddle points (see e.g.B). 
The saddle point equation, Q = (Q + iE)^ 1 , may be solved by going to the basis where Q is diagonal. One finds two 
possible solutions for each of the eigenvalues Ai, . . . , A„ of Q + IE: 
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A ± (£)=±yi-^ + !|. (10) 

This results in 2" distinct diagonal matrices. As Q is obtained from the diagonal matrix through a unitary transfor- 
mation, one finds n + 1 saddle point manifolds, generated by rotations of U (n) / [U (p)U (n — p)] , applied to a diagonal 
matrix with A.; = A_ for i = 1 . r^p and Xj = A+ for j = p + 1 . . . n, where < p < n. Which of these manifolds 
dominates is not obvious a-priorili-a. It is known that the leading, large N, contribution to the DOS is given by the 
simple 'replica symmetric' saddle point (this one is just a point not a manifold, since Q is proportional to 1) with 
p = 00. We show below that other saddle point manifolds are crucial for computations of the 1 /N corrections to the 
DOS, as well as the level correlations. In order to take into account all possible saddle points, one must be careful 
to integrate_Droperly over the manifolds. To achieve this we found it convenient to employ the method introduced 
by T. GuhrE£l in the SUSY context. This method takes advantage of the Itzykson-Zuber integralE3 for the GUE to 
integrate all rotational degrees of freedom exactly, and discusses the saddle points in terms of the eigenvalues of Q. 
A similar method has, been employed by Brezin and Hikami in order to derive and extend the universality of level 
spacing distributionaiS. I— ■ 

By shifting first the Q matrix Q — * Q — iEtJi one obtains: 

ZM{E)=c n (-i) Nn j d[Q]e X p{-jTr(Q-iE) 2 \{detQ) N (11) 

The Hermitian matrix Q is then diagonalized by a unitary transformation, Q = U~ 1 AU with U £ U(n) and the 
diagonal matrix A = diag{Ai, . . . A„}. The volume element transforms like 

d[Q]=A*(A)d[A] dfx(U); 

n 

d[A] = l[dX j , (12) 

where dfj,(U) is the measure of the group U(n) and the Jacobian is given by the square of the Vandermonde determinant 

A„(A)= J] (Xj — Ai) . (13) 

l<i<j<n 



The non-trivial integration flver the group U (n) in Eq. ( pi] ) may be performed using the Itzykson-Zuber integrals 
(see also Appendix 5 of Ref .El) : 

ZW(Ej = Cn(-i) Nn J d[A]A r 2 l (A)(detA) Ar j dfx(U) exp j-y Tr(A - iUEIJ- 1 ) 2 

^c n (^) Nn (^y {n ' 1)/2 ^^Jd[A}A n (A)(dctA) N Cxp L . (14) 



It may seem that the last expression has poles at Ei = Ej, which is not the case. Indeed, the integral over the 
eigenvalues Xj results in a totally antisymmetric function of Ei, which vanishes if Ei = Ej. The reason we introduced 
the diagonal matrix E with all elements different was to regularize properly this fictitious singularity. The next step 
is to take the limit E\, . . .E n —> E. This procedure is described in appendix The result is given by 



Z n (E) = c' n J d[A] A 2 n (A) exp j - N £ A( A, , E) J , 



(15) 



where 



1 AT™ 2 / 2 1 

A(A^) = -(A j - l E f In A, , d n = (-i) Nn - - . (16) 

2 {27r) n j! 



Equation (|T^) may be easily derived already from the first line of Eq. (|l4|) if E is proportional to the unit matrix. The 
reason the longer procedure was presented is to generalize it later for the case of multi-point correlation functions. 
Employing Eq. (^|), one finally obtains for the replicated single-point correlation function 
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s[ n \E)=ic! n /d[A]A^(A) exp ^ —N\]A(Xj,E) 



(17) 



So far all the calculations were exact. One may take now advantage of the large parameter N 3> 1 to perform 
the integrations in Eq. (|l^) by the saddle point method. Differentiating A(Xj, E) over Ay, one finds two saddle point 
solutions for each Xj, which are given by X±(E) defined by Eq. (|l0|). This leads to 2 n distinct saddle points of the 
integral in Eq. ([n]) , each of them may be brought to the form 



A = diag{ A_, . A_ , A + , . A+ } . 

p n—p 



(18) 



There are = n\/[p\(n — p)\] such saddle points for every p, < p < n. Finding dominant saddle points is not totally 
trivial. On the oue. hand, the saddle point action, A± = A(X±(E), E), is such that \A+\ = \A-\ when the energy is 
real (see footnoteta) . Furthermore, one must be careful with the saddle point calculation of the integral in Eq. (|17| ) 
because the preexponential factor, A^(A), vanishes identically at any saddle point (for n > 2). Therefore it should 
be expanded to sufficiently high power of Xj — X± to produce a non-zero result. To this end we introduce variables 
£j describing fluctuations around the saddle point, Eq. (fl8|), defined as: 



(19) 



Xi = X- + d/vN ; i = l...p, 
X 3 = X+ +£j/VN; j=p+l...n 



and the diagonal matrices S_ = diag{£i . . . £ p } and s+ = diag{£ p +i . . . £„}. For any p one may identically rewrite 
A* (A) as 



^(A) = 



p(p-l) + (n-p)(n-p-l) 



II II 

i— 1 j—p-\-l 



Si - 6 



N 



A*(3_)AjU p (5 H 



(20) 



The factor in square brackets on the r.h.s. of this expression is non-vanishing at the saddle point and therefore may 
be substituted by its saddle point value. Expanding the exponent to second order in the deviations from the saddle 
point, one obtains for the replicated correlation function 



S{ n) (E) = ic' n J2 CI [A+ - A_] 2P( "~ P) e -NpA-N{n- P )A + f A _ + (n _ p)x+ _ 



p=Q 



1 



p 2 +(n-p) 2 



/n^A^S^expl-^^e?) ff[ ^A 2 „_ p (H + )exp(-^^ A 
i=i { i=i J J j= P +i [ j= P +i J 



(21) 



where A'± = d^A(X, E)\ \ ± . The two remaining integrals are known as a version of Selberg's integral^, given by 
ff[^ i Al(E^e W i-tJ2&) = m- p2 / 2 n p , with (l p = (2iry/*f[i\. 

i=l { i=l ) i=l 



(22) 



As a result the correlation function takes the form 



s[ n) (E) = (») 



1—nN^-nNA 



n 

+ ^F£[A + -A_] 

p=0 



2p(n-p) 



J\fp(n-p) 



{jA^i)( n -py {JA^y 2 



oP N(A + -A^) 



pX^ + (n — p)X + — inE 



(23) 



where we have introduced the symbol as 



pp = np 



p n—p 

n i! n ?■ 

3 = 1 3=1 
n 

n j! 

3=1 



n r(n _ i+1) 



P^O 



(24) 



4 



and FP =0 — 1. Since the gamma function diverges at any negative integer, F^ >n = 0. Therefore one may extend 
the summation over p in Eq. (^3| ) up to infinity. The resulting expression is suitable for the analytical continuation, 
n — > 0. To continue the F% function one may use the identity T(—z) — —tt/(T(z + 1) sin7rz) with z = j — 1 — n, which 
leads to 



FP = (-l)P(f- 1 )/ 2 



i[mr(j-n); 

3=1 



(25) 



Expanding the sin 7m for small n, one obtains i 7 ^"" — 1; F%z} n = n, whereas F%z? = 0(n p ). As a result only two 
terms with p = and p = 1 survive in the sum in Eq. (23). This is an important conclusion: out of the 2™ possible 
saddle points only one with p = and n with p = 1 contribute in the replica limit. Of course, the argument we have 
given here is somewhat heuristic since the series in Eq. (|2^) is divergent for non-integer n. This should be justified 
more rigorously. Although we have not completed the rigorous proof of this statement, we present some elements of 
the proof in appendix [b|. 

One may easily check that the p = 1 contribution is smaller by a factor l/N with respect to the p = one. We are 
therefore back to the familiar statement that the large N limit of the DOS may be calculated by the replica method 
using a single trivial saddle point for the Q matrix which is proportional to the unit matrix. In particular one obtains 



S 1 (E) {p=0) =i(X + -iE). 

Employing Eqs. (|), (0) one finds the famous "law of semi-circle" for the DOS, vq(E) = V4 — E 2 /(2tt). 
It is instructive to look at the l/N contribution originating from the p = 1 saddle point manifold: 



(26) 



s i( E )(p=D 



1 



1 



,N(A+-A-) 



This leads to an oscillatory correction to the mean DOS of the following form 



(27) 



1 (-l) w+1 1 
SVosc(E) = — \_ E 2 cos[A (29 + sin 29)} = 



1 



N4n 3 v%{E) 



cos 



2irN J is (E)dE + ir 

-2 



(28) 



where sin# = E/2, This is indeed the correct l/N oscillating correction, as can be checked directly using the exact, 
finite N, expression for the DOSEl: 



v{E) 



,-NE 2 /2 




" x ( E\l y) Hn ( E \Iy)~ HN - 1 ( E \Iy) HN+1 ( E \l f 



V2/nN2 N T(N) 

and employing the following integral representation of the Hermite polynomials 



(29) 




N+k+l oo 



dXe -NA(X,E) 



(30) 



Evaluating this integral in the saddle point approximation, one obtains the "semi-circle" as the leading term and 
Eq. (28) as the l/N correction. We do not discuss here the l/N correction to the smooth part of the DOS, which 
may be easily evaluated by expanding near the trivial p = saddle pointESEI. 

Notice that in deriving this result all n+ 1 saddle point manifolds in the Q space were taken into account. However, 
only two of them appear to contribute in the n — > limit. In the leading order in N — * oo only the trivial one Q ~ 1 
remains, whereas the p = 1 is responsible for the oscillatory l/N contribution to the DOS. We shall see below that 
in the case of the two-point correlation function non-trivial saddle point manifolds contribute already to the leading 
order. 

These corrections due to the replica symmetry braking saddle point, which lead to an oscillatory behaviour of the 
density of states in the range A £ [—2,2], also lead to an exponentially small (in 7V)_tail of the density of states outside 
of the interval [—2,2]. This fact was first noticed by Cavagna, Giardina and ParisO, who also showed that the p = 1 
saddle point reproduces the correct exponentially small tail. 
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III. TWO-POINT CORRELATION FUNCTION 



The two-point correlation function is defined as 



S 2 (E, E>) = N-^r(E-H)-^{E>-H)^= N~ 2 dhl ^ , (31) 

where complex energies E and E' have negative and positive imaginary parts correspondingly. Introducing two sets 
of replicas with sizes n and n' respectively to handle each of the logarithms in this expression, one obtains 

S 2 (E,E>) = lim ^-Si n+n '\E,E>); S< nW) = JLJ?L^Z{E)»Z(E>y . (32) 

n,n'^o nn' N z oEoE' 

We introduce again, for the sake of regularization, the function Z^ n+n ' (E), where E now, and in the rest of this section, 
is a diagonal (n + n') x (n + n') matrix of the form E — diag{i?i, . . . E n , E n+ \, . . . E n+n r}; the limit E±, . . . E n — > E; 
E n+ i, . . .E n+n r — > E' will be taken at the appropriate stage. The next steps are exactly identical to those of the 
previous section, up to the change of n into n + n': averaging over H , decoupling with the (n + n') x (n + n') matrix 
field Q and integrating over the group U(n + n') using the Itzykson-Zuber integral. This leads to the following result 
for the function Z {n+n '^E) (equivalent to that of Eq. @): 



n+n' 



where A is a diagonal (n + n') x (n + n') matrix containing the eigenvalues of the Q matrix. Once again, since the 
integral is a totally antisymmetric function of E^ there are no poles at Ei = Ej. The next step is to take the limit 
Ej — > E for j = 1 ... 7i and Ey — * E' for j' = n + 1 . . . n + n' . The corresponding limit is calculated in the appendix 
[A[, the result is: 

, / 

I! II (V-AjO f „ n+n' ) 

Z^E)Z-'{E')=d n d n , yd[A]A^(AW)A^(A("0) 3 -j^+ 1 j;)]wn) expj-JV^^Aj.^-^E A ( X f> E V > ( 34 ) 

where AW = diag{Ai . . . A„} and A^™ ) = diag{A„ + i . . . A„ + „'}. Differentiation with respect to E and E' gives the 
replicated two-point correlation function. Non-trivial correlations exist only in the range where E — E' is of the order 
of the level spacing, namely 1/N. We introduce thus scaling variables e = NE and e' = NE' and restrict ourselves to 
the vicinity of the center of the band |e|, |e'| -C N. The large N limit of the correlation function at fixed e, e' is given 
by: 

n n+n' 

, n n (a, 

S { r n '\e,e') = c ;4,iV«™7 d [A]A ? 2 l (A(«))A^(A(™'))^g± 1 



[i(e'-e)\ 




n n+n' 
j=l j'=n+l 

In the large N limit the correlation function may be evaluated using the saddle point approximation. The saddle 
points of the integral in Eq. ( |35| ) are given by \j = ±1, for j = 1, . . .n + n' , which is the zero energy limit of 
\±(E) discussed in the previous section. Altogether there are 2 n+n distinct saddle points. Each one of them may be 
parametrized in the following manner 

A = diag{-l, . . . - 1, +1, . . . + 1, +1, . . . + 1, -1, . . . - 1} , (36) 




where < p < n and < p' < n' . Given p and p', there are C^C^, such saddle points. As before, one must be careful 
with the saddle point calculation because of the vanishing prefactors. The method we shall follow is the same as in 
the previous section. One introduces variables describing fluctuations around the saddle point 



G 



Ai= -1+&/VN; i = l...p, 

Xj = + 1+£j/\/N; j = p + l...n, 

Xi' = + l+Zl'/VN ; i' = n+l...n + p' , 

Xj, = - ; j' = n + p' + 1 . . . n + ri . 



(37) 



and groups them into diagonal matrices ^_ = diag{£i . . .£ p }, S + = diag{£ p+ i . . .£„}, E' + = diag{£^ +1 . . .^ n+P '}, 
S'_ = diag{£^ + , +1 . . . Cn+n'}- F° r l ar g e ^ the determinants are decomposed as: 



p(p-l)/2+0-p)(n-p-l)/2 

2 f("-f)A p (S_)A n _ p (H + ), 

p'(p'-l)/2+(n'-p')("'-p'-l)/2 



(-2f'(»'-^)A y (HV)A„,_ y (H'_). 



(38) 



Finally the remaining factor in Eq. ([55|) takes the form 

2pp (_2)(™~p)("' - p') 



n n+n 



n n (x r - a 7 ) _ ._ 

11 11 31 (^/N)p(n'- P ') + (n-p)p' 

j=lj'=n+l vv I 



p n+n 



ii n 

1 = 1 j'=Tl+p' + l 



n n+p 



n n 



(39) 



Expanding the exponent in Eq. (Bq) around the saddle point, Eq. (139), 



iV 



iV 



4X 



2VX 



/ 1 



\ N 3 /< 



(40) 



one finds that the integrals over £j and £j may be expressed as I p . n '-p> u(e — e')/2\/]V^ x In-p,p' — e')/%V~Nj , 
where the function J riS (a) is a generalization of Selberg's integral, defined as: 



7 r , s (a)= / d[X]d[y]A r (X)A s (y)A r+a (A:©y)expJ -^fe-a) 2 -^^ I 
^ [ j'=i fc=i J 



(41) 



Here X, Y and X © Y are diagonal matrices: X = diag{xi . . . x r }, Y = diagjyi . . . y s }, and X © Y 
diagjxi . . . x r yi . . . y s }. We show in appendix ^| that this integral is given by 



(42) 



where il r is the usual Selberg integral, Eq. (g2|). Up to an overall constant factor which goes to one in the limit 
n, n' — > 0, one has thus: 



r<pr<P j 



i(e-e') 



2VN 



i(e-e>) 



2VN 



n' 



i(e'-e) 



N 



p(n' —p') + (n—p)p 



Grouping all the terms, one finally obtains 

(— i)pp' 2~ 3p2 ~ 3 p' 2 '+ i pp' j\j( n ~p+p')( n ' '-p' '+p) 



p=0p'=0 



Ul e > _ e Xlnn' -p(n' -p')-(n-p)p' 

inn' 



2~(p-p') 2 



,2i(p e —pe)-\-i-!rN(p—p ) 



(43) 



n - 2p n' - 2p' - ^ r (n - 2p + n' - 2p') ± ^-L 



(44) 



where we have omitted an inessential factor const°^ n \ 

Employing the fact that F£ >n = F%, >n = 0, one may extend summations over p and p 1 to infinity and then perform 
the analytical continuation, n,n' — > 0. Due to the properties of the F-symbol (cf. Eq. (p5|)) only the terms with 
p = 0, 1 and p' — 0, 1 contribute in the replica limit. We need to evaluate the contributions of these four saddle 
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points. Recalling that F® = F®, = 1, one finds for the p = p' = contribution to the two-point correlation function, 
S 2 (e,e')=lim n , n ^o(nn')- 1 Si n+n ' ) : 



S 2 (u,) (p=pl=0) =[l-^) , (45) 

where u> = e — e'. This is the result obtained in the perturbation theory around the usual saddle pointfl There is, 
however, an other contribution to S 2 originating from the saddle points with p = p' = 1. It is easily computed, using 
the fact that F ^ = n and Fh = n', and is equal to: 

&M (p=p , =1) = (46) 

It is easy to check that the saddle points with p = 0, p' = 1 and p = 1, p' = lead to 1/N oscillatory correction 
to the disconnected part of the two-point correlation function, cf. Eq. (pgj). Adding together the two leading terms, 
Eqs. (pf5|), ((46|), one obtains the final result for the two-point correlation function : 

1 — e~ 2iu> 

S 2 {io) = 1 - 2 = 1 - 2iLU- 2 e - luJ sinio. (47) 



Although this is the exact resultQ'llHl9'E3 for the GUE for any u <C N , the way it was derived here justifies it only 
for 1 <C to <C N. This is because of the terms which were omitted in the expansions (H^), (|39|). Although they look 
superficially as being of order 1/^/N, they can, in fact, contribute. A careful examination shows that the generalized 
Selberg's integral is just the leading large lo (to ^> 1) contribution. We have thus computed only the leading term 
at large 10 for each saddle point, but corrections could be incorporated systematically. Accidentally the obtained 
expression appears to be exact down to zero to. A similar situation was already encountered in a SUSY approach, 
when_the soft modes integrals were calculated with the saddle point method including an additional non-trivial saddle 
pointEJ. In this sense our new saddle point with p = p' = 1 is a close analog of the SUSY saddle point of Andreev 
and Altshulery. 



IV. DISCUSSION OF THE METHOD 



Let us add some comments on the relation between our calculations and the saddle point evaluation of the sigma 
model, Eq ^ with E = diag{_El n , E'l n '}. Looking for the saddle point solution in the form Q = C/ _1 A[7, one finds 
the solution in the form 

o-(va)»(:;). w 

with A being a diagonal matrix, obeying A 2 + iEA —1 = 0, and arbitrary V E U(n); V' E U(n'). For a diagonal 
matrix A of the structure given by Eq. (|36"|), there is a set of rotations belonging to the coset space 

m x -W (49) 

which leave the action invariant, while changing the saddle point matrix, Q. As a result, there is a continuous saddle 
point manifold, which contains true zero modes of the functional integral, Eq (Js[) . In addition there are usual "soft 
modes" with masses of order |e — e'\ <C N. Notice, that there are no zero modes around the trivial saddle point 
p = p' = 0. The saddle point, Eq. (p6|), contains n — p + p' components which are +1 and n! — p' + p components 
which arc —1. Therefore, out of the total (n + n') 2 fluctuation directions (n + n') 2 — 2(n — p + p')(n' — p' + p) 
are massive with the mass JV, whereas the remaining 2(n — p + p')(n' — p' + p) degrees of freedom are split between 
[n 2 -p 2 -(n-p) 2 ] + [n' 2 -p' 2 -(n'~p') 2 ] = 2{j>(n-p)+p' (n' -p')) zero modes (cf. @) and 2nn'-2p(n'-p')-2(n-p)p' 
soft modes with the mass |e — <C N. The integrals over the zero modes must be calculated exactly giving rise to 
the volume of the coset space (fly). In the regime 1 <C |e — e'|, the integrals over both massive and soft modes may 
be evaluated in the Gaussian approximation giving rise to factors iV -1 / 2 and |e — e'| -1 / 2 in the number of modes 
power. This is precisely the structure of Eq. (|44|), where the factor FPF%, is proportional to the volume of the coset 
space (fl9]). The advantage of our method is an easy control over combinatorial factors, coefficipHts etc., otherwise it 
is equivalent to the Gaussian evaluation of the functional integral, Eq (Q), similar to that of Reffil for the SUSY case. 
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Our method can be easily generalized for the higher order correlation functions. For example, calculations of the 
three-point function S^(E,E' ,E") with, say, E having negative and E',E" positive imaginary parts, lead to the 
triple sum over p, p', p" analogous to that of Eq. ([44|). Again only the terms with p,p',p" — 0, 1 contribute in the 
replica limit. One may easily check that the correct result for the connected GUE three-point correlation function 
follows from the p = p' = 1; p" = and p = p" = 1; p' = terms, whereas all other possible combinations, including 
p = p' = p" = 1, appears to be small in powers of I/TV. 

V. CONCLUSION AND PERSPECTIVES 

There are several questions raised by our computation. An obvious one is to have a more rigorous derivation of the 
analytic continuation of the g(x, n) function at small n (see appendix ^). Extending this approach to bosonic replicas 
and to other random matrix ensembles are among other open problems. 

We have presented here what we believe is the first consistent application of the replica method to the RMT. Our 
computation reconciles the fermionic replica result with the previous approaches. The point of this paper is not to 
challenge these previous approaches. The results which we have derived here have been weiLknown for years, and in 
fact there exist in the litterature much stronger results on level spacing universality (see e.g.Ej and references therein). 
The a model representation itself has proven very successfull when used with the SUSY method: in problems of 
random energy levels, the SUSY technique has been very well [developed and has allowed to derive many results 
in various problems of solid state and nuclear physics (see RefsEd'EI). We think that our result has two interesting 
aspects: the mathematical consistency on the one hand, and the possibility to use these ideas for a study of disordered 
interacting electrons. 

The previous situation in which the replica approach was considered as ill was not satisfactory from the mathematical 
point of view. Furthermore the replica method is known to be highly successful in other probleaas such as the statistical 
mechanics of classical disordered systems (see Ref.E3 for a review) and localization theoryE3, and its failure in the 
simple problem of eigenvalue correlations seemed strange. In this respect we would like to comment about the replica 
symmetry breaking (RSB) which we have found. While the cr-model formulation seems to involve a n x n matrix 
order parameter, similar to the one which has been discussed for instance in spin glass problems, the symmetry groups 
are very different. In spin glasses the symmetry group of the replicated system is just the permutation group, while 
in the cr-model there ia a larger symmetry group: in our case some version of the unitary group (depending on the 
type of correlation one computes, and whether one would use commuting or anti-commuting replicas). Integrating 
over the angular variables has left us with an order parameter (the set of n eigenvalues of the Q matrix) which is a 
vector in the replica space. Therefore the pattern of the replica symmetry breaking which we have found is much 
more reminiscent of the 'vector RSB', discussed in the study of random-field-like problems at low temperaturetj, 
rather than the hierarchical RSB scheme, which describes the spin glass mean -field theory. The vector RSB may 
be traced back to the existence of several distinct ground state configurations in a problem. In all cases studied 
so far, there is an infinity of RSB saddle points, which contribute to the partition function. At the same time the 
SUSY approach cannot address these problems, because it is unable to estimate a sum over ground states, but rather 
computes a topological invariant, given by the sum over all saddle points, weighted by the parity of the number of 
unstable directionsEJ. In the RMT the situation is much simpler: SUSY is exact, and there is only a finite number of 
saddle points contributing in the vector RSB (two saddle points in the DOS computation). These facts are certainly 
related, and it would be highly desirable to understand better their connection. 

The SUSY method relies crucially on the fact that the original action (as in Eq. (|5])) is quadratic in the field 
variables. In the application to electronic system, it is thus restricted to non-interacting electrons. The renlica 
method does not have such a limitation, and it is therefore capable to address problems of interacting electronsc3. It 
would be very interesting to see whether the new saddle points which we have found have some implications in the 
theory of interacting electrons in disordered media. 
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APPENDIX A: 



We first evaluate the integral given in Eq. (|l4|), in the limit Ej — > E, j 6 {1 . ..n}. Denoting Ej = E + r/j, one 
rewrites the integral in the following form 



r [ n n N - 1 

C(Ei ,...,E n )= d[A]A n (A) exp I —N ]T A(Aj , E) + nJ^ *j (%) - y £ ^ f 

[ i=i i=i j=i J 



(Al) 



where A(\j, E) is defined by Eq. (|16|). Expanding the term exp{iV Aj-(irjj)} in series, one obtains: 



. AT oo iyfci + ...fe n 

C(i?i,...,K) = exp<j- y E^r E fc I k j (iVi) kl ---(iVn) k "T kl ... K , (A2) 

j I fci,...,k„=i 1 " " '"' 



where the tensor T is a function of E defined as: 



T, 



:i ... fc „ = J d[A]A„(A)A^ 1 ...A^exp j-A^A(A,,£) j . (A3) 



Since the Vandermonde determinant is antisymmetric in A.,-, the tensor T is fully antisymmetric: for any permutation 
7r of the n indices, one has 

^...k, = Ji, w .j, w Si , (A4) 

where SV = ±1 is the signature of the permutation. In particular, one notices that X^.../^ vanishes whenever two 
exponents ki and kj (with i ^ j) are equ al. We are interested in the leading behavior of £ when all r\j go to zero 
simultaneously. From the expression ( |A2[) and the antisymmetry of T, it is clear that the leading term is of order 
^n(n-i)/2 anc j j g obtained whenever k\ = 0, k% = 1, . . . , fc n = n — 1, or any permutation of the integers 0, 1, ... n — 1, 
and the exponential prefactor may be neglected. The leading order in 77 may be written as a sum over all permutations 
7r of the ensemble {0, . . . , n — 1}: 

C(£a, . . . , E n ) ~ T 01 ... n _! aK«-i)/2_L_ £ 5,(^)^(0) . . . (^("-U (A5) 

11.7=0 7T 

In the sum over permutations one recognizes the Vandermonde determinant of the ir)j, which is equal to A n (iE) 
and thus cancels the corresponding factor in the denominator of Eq. (p"4|). As for the value of 2oi... n _i, it may be 
rewritten, using again the antisymmetry of A„(A), as: 

Toi...„-i = ^ £ S - J <A] A„(A)A? (0) . . . A^"- 1 ) exp ^~ N J2 A ^ > E ) 

= ±j rf[A]A^(A)exp J -NJ2M^,E) 1 . (A6) 



Therefore, in the limit rjj — > 0, C behaves as: 



C(E 1 ,...,E n )^N n ^- 1 ^ 2 ^— -A n {iE) J d[A]A2(A)exp|-7VEA(A J ,i?)| . (A7) 

This establishes expression (|l5|). 

We now evaluate the n + n dimensional integral appearing in Eq. (|33|), in the limit where Ej — ► E for j = 1 . . . n 
and Eji — > E' for j' = n + 1 . . . n + n'. The procedure is exactly the same as was explained above for the one-point 
function. We shall sketch it briefly. One writes Ej = E + r/j and Eji = E' + rj',,_ n in terms of which the integral reads: 
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. in n+n' I 

C(E 1 ,...,E n+n ,)= d[A]A n+n ,(A)expl-Nj2 A (^,E)-N £ A(X f ,E')\ 

J [ j=l j'=n+l J 



in n N N \ 

j=l j'=l J= l j' = l J 

We expand the terms exp{iV^ ■ Xj(irjj)} and exp{iV^Z ■, A, i+ y (z^., )} in series, and notice that in the limit where all 
the rjj and rj-, go to zero simultaneously and independently, the leading contribution, of order (f7)"( ra_1 )/ 2 (r/)" ( n _1 )/ 2 
is obtained when powers of irjj (resp. irj'j) span the ensemble {0, . . . , n — 1} (resp. {0, . . . ,n' — 1}). The result can 
be written as a sum over all permutations tt of the ensemble {0, . . . , n — 1}, and all permutations tt' of the ensemble 
{0, . . . , n' — 1}, in the following form: 

(n-1 n'-l \ 1 
J=Q j'—0 J 7T7T' 

(A9) 

where u is equal to: 



d[A]A n+n ,(k) A? A* . . . Ar 1 A° +1 A^ +2 . . . a££ cxp I —N ^j4(Aj,E) - JV £ A(A y , £') I 

[ 3=1 j'=n+l J 



(A10) 



Exactly as above, this integral may be rewritten by permuting separately the dummy integration variables Ai, . . . , A„ 
and An+i, . . . , A n _|_„/. One thus obtains in the limit where all r]j, tj'j, — ► 0, in the notations of Eq. (p4|): 



n n 



/ d[A] A n+n ,(A)A n (A( n ))A„,(A( n '))exp^ -iV ^ A(A„ £) - iV £ A(A,v,£')L (All) 

[ 3=1 j'=n+l J 



This establishes expression (B4f). 



APPENDIX B: 



To prove the statement that only the terms with p = and p — 1 contribute to the analytic continuation of 
Eq. fl25| ) at small n, let us study the function g(x,n) — J2p =0 F%x p . Using the Gaussian decomposition of the 

factors containing exponents of p(n — p) and p 2 , one may show that the correlation function sl^j is deduced from the 
knowledge of g(x,ri), where a; is a complex number, with a modulus slightly smaller than onee3. Using the fact that 
F% +1 = FP r(p+ l)/r(n — p), one finds that for integer n the function g(x, n) satisfies the following integral equation: 

oo 

f f du e u 

g(x, n) = 1 + x I dr e T — —g(xTu,n), (Bl) 
o c 

where C is the contour in the complex u plane used for definition of the function 1/T(z): it goes around the negative 
real half axis, starting from — oo to with a small positive imaginary part, turning around and getting back to — oo 
with a small negative imaginary part, it thus passes around the cut of the l/u n function for n non-integer. This 
integral may be probably used to define the function g(x,n) for an arbitrary n, although some further study of this 
statement is needed. Here we are interested in the behavior of g{x,n) at small n. Writing the first two terms in the 
small n expansion as g(x, n) = go(x) + ngx(x) + . . ., one finds that go and g\ satisfy the following equations: 
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oc 

g (x) = l + x J dre~ T J -^-e u g {xTu) 
o c 

oo 

gi(x) — x J dre~ T J — e u [gi(xru) — g (xru) Inu] . (B2) 
o c 

Assuming that go and g\ are analytic in a certain domain T> near the origin, one can compute them inside this domain 
by series expansion in powers of x. This leads immediately to go — 1 and g\ — x, which gives exactly the same answer 
as our heuristic arguments given in the text. To complete the proof one has to find out a shape of T>. We believe that 
T> is the part of the complex plane restricted by the unit circle, but we have not been able to prove it 

APPENDIX C: 



In this appendix we prove that the generalized Selberg integral, Eq. (|4l|), is given by the expression (42). The proof 
consist of the two steps: 

(i) We shall prove that the series expansion of the integral in powers of a starts as a m with m > rs. To this end we 
rewrite the integral as 



IrA a )= I d[X]d[Y] A r (X)A s (Y)A r+s (X ®Y)exp { -^a;| -^y 2 k + 2a^x J - ra 2 } (CI) 

3=1 k=i j=i 

and notice that the integrand is equal to A r (X) times a totally antisymmetric function of the Xj. This allows one 
to substitute in the integrand A r (X) by rl x\x\ . . ,x* . A similar observation for the y variables allows one to 
substitute in the integrand A S (Y) by si 2/i 3/2 • ■ ■ 2/J 1 giving: 

I r , s (a) = r\s\e~™- f d[X]d[Y] A r+s (X © Y) \x\x\ . . . x^ 1 y\y\ . . . y^ 1 } exp J - ]T x) - £y» + 2aj^ x 3 \ (C2) 

J [ 3=1 fc=l 3=1 J 

The integrand is the product of a term which is totally antisymmetric in all the r + s integration variables times the 
factor [x\x\ . . . x r r ~ x y\y\ ■ ■ ■ yl, 1 } exp{2a^ J - Xj}. In this factor one can expand the exponential in a power series in 
a. Whenever there are two of the r + s variables appearing with the same power, the the integral is zero as can be 
seen by permuting these two variables. The first non-zero contribution appears, thus, when the power series generates 
a power like y\y\ . . . yl" 1 x^^ x^ s+1 ^ . . . ir' s+r , where 7r is any permutation of the integers s, s + 1, . . . , s + r — 1. 
Such terms appear when the exponential is expanded to the order a rs . This shows that the series expansion of the 
integral in powers of a starts at least with the order a rs . 

(ii) We demonstrate now that I TtS is a polynomial in a of degree less or equal to rs. Shifting each Xj to Xj + a and 
splitting the factor A r+S (A © Y), one may rewrite the integral as 

/[ r sirs 

d[X]d{Y] A 2 r (X)A 2 s (Y) exp <^ - £ x) - ]T y\ \ l[(y k - x, - a) . (C3) 

[ 3=1 fc=l J i=lfc=l 

This shows that the integral is a polynomial in a of a degree less or equal to rs. 



As a result of (i) and (ii), I r , s { a ) mus t be proportional to a rs . Using Eq. (C3), one finds that it may be expressed 
as a product of two usual Selberg integrals, Eq. (|2|), as: 

IrAa) = (-a) rs I d[X]d[Y] A 2 r (X)A 2 s (Y) exp J - £ x 2 - £ y 2 k i = (-a)-2-^ 2 + s2 )/ 2 ^ r ^ . (C4) 

3=1 fc=l I 
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